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Abstract 

We present the simple and direct proof of the determinantal formula for the scalar 
O ■ product of Bethe eigenstate with an arbitrary dual state. We briefly review the direct 

calculation of the general scalar product with the help of the factorizing operator and 
the construction of the factorizing operator itself. We also comment on the previous 
Q ! determination of the scalar product of Bethe eigenstate with an arbitrary dual state. 
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1. Introduction. 



> 

CN \ One of the most important open problems in the theory of quantum integrable models is 

the calculation of the correlation functions. In the framework of the Algebraic Bethe Ansatz 
method |I[ the problem is the combinatorial complexity of calculations due to the structure 
of Bethe eigenstates. In spite of these difficulties many important results have been obtained. 



Among them one can mention the derivation by Korepin the Gaudin formula for the norm 
of Bethe eigenstates @ and the calculation of the scalar product of Bethe eigenstate with an 
arbitrary dual state ||, which leads straightforwardly to the determinant representation of 
the formfactors of local operators (see j| for the case of spin- 1/2 chains and 0, || for the 
case of Bose-gas with S - function interaction). Using the Algebraic Bethe Ansatz and the 
solution of the quantum inverse scattering problem the authors of ref . J7| obtained the multiple 
integral representation for the correlation functions found previously using the other methods 
(see for example || and references therein). The concept of factorizing F - matrix introduced 
recently by Maillet and Sanches de Santos || following the concept of Drinfeld's twists in his 
theory of Quantum Groups and construction of the Monodromy matrix in the F -basis ||, 0, 
||T0| allows in some respect to simplify the calculations in comparison with the general theory 
of scalar products developed previously 0, fllT |. 



At the same time only the very limited number of physical results for the correlation 
functions (or formfactors) have been obtained from the first principles. So the developing of 
new calculational methods and better understanding of the underlying mathematical structure 
for the expressions for the correlators is required. For example it is still not completely clear 
if it is possible to obtain the determinantal expressions for the correlators without using the 



auxiliary dual fields [fj , [1 1 



The aim of the present letter is to present the simple and direct proof for the expression of 
the scalar product of Bethe eigenstate with an arbitrary dual state. First, to fix the notations 
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and derive the expression for an arbitrary scalar product, we briefly review the results of 
ref.jllj on the construction of the factorizing operator (Section 2), and simple calculation of 
the general scalar product (Section 3). Next, we present the new version of the proof given 
by Slavnov || for the scalar with Bethe eigenstate. Finally in Section 4 we present the direct 
and a simple proof of this formula. Our results could shed some light on the possibility to 
obtain the determinantal expressions for various correlators. 

2. Construction of the factorizing operator. 

We consider in this letter the XXX or XXZ spin- 1/2 chains of finite length N. In the 
present section we diagonalize the operator A and introduce the factorizing operator. Let us 
fix the notations: the normalization of basic S - matrix, the definition of monodromy matrix 
and write down the Bethe Ansatz equations. For the rational case (XXX- chain) the S- matrix 
has the form S^^i,^) = ti — h + f]Pi2, where P±2 is the permutation operator. In general 
(XXZ) case it can be written as 

/ a(t) \ 

c{t) b(t) 

b{t) c(t) 

V a(t) / (1>i 



Sl2(£l, £2 



t — ti — to. 



One can choose the normalization a(t) = 1 so that the functions b(t) and c(t) become b(t) = 
4>(rj) I (f>(t + 77), c(t) = (p{t) /(pit + 77), where 0(t) = t for the isotropic (XXX) chain and 0(t) = 
sinh(t) for the XXZ- chain. With this normalization the S'-matrix satisfies the unitarity 
condition 5 , i 2 (ti, t2)5 , 2i(^2, ^1) = 1- The monodromy matrix is defined as 

T (t, {£}) = S 10 (£ 1 ,t)S 20 (&,t)...S m ({; N ,t), 

where £j are the inhomogeneity parameters. We define the operator entries in the auxiliary 
space (0) as follows: 



(8\T\a)-( A{t) B{t) 
mTola) -{c(t) D(t) 



I ; a,/3 = (l,2) = (TU). 

a/3 



We denote throughout the paper (t; I) = (1; 0) so that the pseudovacuum (quantum reference 
state) |0) = |{00...0}at). The triangle relation (Yang-Baxter equation) reads: 

5'l25'l3>S'23 — >S'235'l3>S'l2 ) Rq ,TqTqi = TqiTqRqqi] Rqqi = Poq' Sqq' ■ 

The action of the operators on the pseudovacuum is: A(t)\0) = a(t)\0) (a(t) = l\ a c(£ Q — t)), 
D{t)\0) = |0), C(i)|0) = 0. The Bethe Ansatz equaions for the eigenstate of the Hamiltonian 
Yl^Li B(ti) |0) and the corresponding eigenvalue of the transfer - matrix Z(t) = A{t) + D(t) are 

M M 

<ti) = n - uwiti - t a ))-\ A(t, {t a }) = a (t) n c~\t a - 1) + n - 
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where t a are the solution of Bethe Ansatz equations. 

In order to construct the operator O = Ox...n which diagonalizes the operator A(t) 
{0~ l AO = diag(A)), let us first construct the eigenfunctions of the operator A. One can 
do it in two different ways. First, note that A - is a triangular matrix in a sense that it 
makes particles (the spin-up - coordinates) move to the right on the lattice 1...N. Thus its 
eigenvalues are coincide with its diagonal matrix elements and therefore are charactarized by 
the set of integers (rii, ...Um) - the spin-up positions. Let us denote this eigenfunctions by 
\<p(ni, ...Um))- Clearly, |{00..0}tv_a/{H..1}a/) is an eigenstate of A{t). Therefore, considering 
the permutation 

Sio-'-SWo S\Q...SNoS nM Q...S ni Q {rii < n 2 < ... < um), 

we realize that 

\(f>(nx,...n M )) = T ni T n2 ...T nM \ni,n 2 , ...n M ) (1) 

where 

is an eigenstate of the operator A(t). Note that if we modify the given permutation inter- 
changing rii n j the operator in the right-hand side of (|T|) modifies but the state \4>({n})) 
remains the same. For example, for M = 2: T n2 T' ni = T ni T n2 S~^ ni , where n\ < n 2 and the 
prime means the absence of the term S nin2 in T ni . Since S nin2 \nx,n2) = \ni,n 2 ) the state 
remains the same. 

The second way - is to consider the state 

|0(n l5 ...n M )) = B(£ ni )B{£^)...B(£ nu )\ti), (m ? n ). (2) 

Using the fundamental commutation relation: 

A(t)B(q) = -±-^B{q)A{t) - b ^^B(t)A(q) } (3) 

and the fact that for the pseudovacuum state A(£j)|0) = 0, we find again that |0({w})) 
{{n} = {nx, ...n M }) is an eigenstate of A(t) with the eigenvalue A^ n ^ n y(t) = Y[ a ^n k — t). 

One can see that the states \4>{{n})) ([[]) and (0) are coincide. This can be seen using the 
identity 

B(U) = T n (0\S w ...S n - lfi Bao\l) (4) 

and ordering the (commuting) operators B(£ ni ) in eq.(|2|) according to the prescription nx < 
n 2 < ... < riM, so that the second operator in the last formula simply creates the particle 
(spin- up) at the site n with the amplitude equal to unity. 

Now we can introduce the diagonalizing operator O = Ox...n ( we show later that it is also 
the factorizing operator). Let us define the operator O such that 

\H{n})) = d({n})\{n}) = d\{n}), 
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where 0({n}) is given by the equation (p. Clearly the operator A F (t) = 0^ l A{t)0 - is 
diagonal. Now we immediately get the simple formula for the matrix elements of the factorizing 
operator: 

{m}{n} = ({m}\B(UB(U)-B(^n M )\0): (5) 

where {m} = {mi, ...tyim}- From this expression (see (|J)) one see that O is the triangular 
matrix (in the same sense as the operator A((t)). Let us show that the operator O is invertible 
and construct the inverse operator. To find this operator we define the operator O analogously 
to the previous case as (<p({n}) \ = {{n}\0 where the corresponding dual states and the matrix 
elements of the operator O are: 

...n M )\ = (0\C{U)-C(Lm), {m}{n} = (0|C(U)C(U)...C(^J|{n}). (6) 
Let us calculate the following scalar product 

(0({m})|0({n») = ({m}\dd\{n}) = (0\C(U)~-C(U„)B(U)--Mtn„)\0). (7) 

This scalar product can be calculated using another well known relation for the commutator 
[B(q),C(t)] following from the Yang-Baxter equation, and using again that for any site i, 
A(£j)|0) = 0. Moving the operators A and D to the right, repeating consequently the relation 
(0) and the analogous relation for D (which differs from eq. (|^) by the interchange t <-> q in the 
coefficients between the products of the operators) and using the equations C(£ n )B(£ n )\Q) = 

Ha^n c(Ca—Cn) |0), we obtain that the matrix / = 00 is diagonal. The corresponding diagonal 
matrix elements can be found either using the procedure mentioned above or, which is the 
simplest way, using the representation @ (and the same for the operator C(t)): 

f(n u ...n M ) = Y[( J] c(e«-U))- ( 8 ) 

Thus we obtained the inverse matrix O^ 1 : 

do = /, 6- 1 = r x b. 

Before proceeding with evaluation of the matrix elements of the other operators in the new 
basis, let us mention some usefull properties of the operator O, and prove that, in fact, it is 
the factorizing operator in a sense of the definition 0. First, O and O -1 are the triangular 
matrices (upper triangular as A(t)). Second, the pseudovacuum state is an eigenstate of O 
(O -1 ) with the eigenvalue equal to unity. In general we have the following equations for 
arbitrary number of particles n: O|{00..0}jv_n{ll..l} n ) = |{00..0}jv_ n {ll..l} n ), the similar 
formulas for the operator O and the same formulas for the inverse operators. From this 
formulas one can already suspect that O is the factorizing operator. Indeed for the particular 
permutation (fl|) cr({n}) the factorizing condition is represented as 0(O tr " n ")~ 1 = 0({n}), 
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where 0({n}) = T ni ..T nM and is fulfilled at least for the state |{n}) since due to the last 
formulas O ff « n ») _1 |{n}) = |{n}). The rigorous proof goes as follows. We construct the 
operator that acting on the state \{n}) produces the state 0({n})\{n}) . It is easy to see that 
the operator that fulfills the above requirement is: 

12 ... N = F 1 F 2 ...F N , F i = {l-n i )+T i fi i} (9) 

where hi is the operator of the number of particles (spin up) at the given site i. The operators 
Fi entering @ do not commute and their ordering in eq.@ is important. To prove the 
factorizing property of this operator it is sufficient to consider only one particular permutation, 
say the the permutation (i,i + 1) since all the others can be obtained as a superposition of 
these ones for different i. One can show ||, |TD] that 



6 = s i+1>i d^ i+1 \ (io) 

Evidently, for any transmutation o G SV we will obtain only one operator R a x N (the operator 
constructed from S- matrices) on the left of the operator O. Thus O is the factorizing operator 
in a sense of ||. 

Let us calculate the matrix elements of B(t) and C(t) - operators in the F - basis: 
B F (t) = FB^F^ 1 = _1 .B(£)0 (and the same for C(t)). The general scheme to perform 
the calculations is to use the formalism developed previously, which leads to the following 
equation: 

B(t)B(U)--MUM = T,( B (^) B (U)---B(Uj\0))<f>(x,t,{n}), (11) 

X 

where 4>{x,t, {n})- is exactly the matrix element of the operator B(t) in the new basis. To 



get the single term in the sum in eq.(|Tl|), we act by the operator A(£ x ) (x ^ njt) at both sides 
of this equation. Using again the property A(£j)|0) = and eq.@ for the left-hand side of 
(|TT| ) we get the single term with B(£ x ) for the right-hand side, which can be evaluated using 
eq.@ and the formula A(£ X )B(£ X )\Q) = Y\ a -L x c{^ a — £ X )B(£ X )\0) , which can be proved by 
direct computations. After the cancellation of similar terms at both sides of eq.([Ll]) we get 
the matrix element 4>(x,t, {n}) and finally obtain in the operator form: 

B F (t) = 5>j m x - 1) n f ~ c{ia " t)m ° " ^ )rl ' 01 1 nk ). (12) 

x a^x V 1> Ol — Uk J 

With this expression one can prove the equation fli B F ((, ni )\0) = \{n}) which is consistent 
with the formulas of the previous section. For the operator C(t) proceeding in a similar way 
and using the relation for the commutator [B; C], we get after some algebra 
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The operators flT2|) and ( |13"D are quasilocal i.e. they describe the flipping of the spin on a 
single site with the amplitude depending on the positions of spins on all the other sites of 
the chain. The operator D F (t) can be found using either the same method or the quantum 
determinant relation and has a (quasi)bilocal form. 

3. Calculation of the scalar products. 

In this section we use the developed formalism to obtain the expressions for the correlation 
functions for the spin chains. We use the factorizing operator and, in particular, the expression 
for its matrix elements (|5|) to obtain the expression for the general correlation function (scalar 
product) 

Sm({A},{*}) = (0\C(X 1 )C(X 2 )...C(X M )B(t 1 )B(t 2 )...B(t M )\0), 

where {A} and {t} are two arbitrary sets of parameters (not necessarily satisfying the Bethe 
Ansatz equations). The correlation function can be represented in the following form 

(0|Tr 0l ,...o 2M (a ( J;...a+ / (7 - /+i ...(7 02M To 1 . . .T 0aM ) |0) 

The auxiliary spaces Oi, . . . 2 m can be considered as a lattice consisting of 2M sites with the 
corresponding spectral parameters Ai, . . . Am, t±, , . . .tu- Rearranging the basic S- matrices 
entering the product of the monodromy matrices we arrive at the operator 

Ti . . . T/v, 

where the new monodromy matrices act in the auxialiary space instead of the original quantum 
space: 

T n = S n o 1 S n Q 2 . . . S n o 2riI . 

Obviously using this matrices the correlator can be represented as the following matrix element 
in the new quantum space Oi, ...0 2 m : 

({00..0} A/ {ll..l} M |i(6) • • • i(6v)|{ii..i} M {oo..o} M ) 

(we use the symmetry A(t) D(t), <->• 1 here). Transforming the operators to the F 

- basis we find at £j = 0: 

Sm = E({00..0} m {11..1}m|O|M)(MI (A f {0))) N |{n})({n}|O- 1 |{ll..l} M {00..0} M ). 

M 

The sum is over the states labeled by the positions of M particles {n} = n\...UM on a 
lattice consisting of 2M sites with the inhomogeneity parameters Ai . . . Xm , t± . . .t M . We use 
O^ 1 = f~ l O and the representations (|5D, (|j) for the matrix elements in the last formula. Let 
us denote by /xi, . . .\iu the parameters (from the set ({A}, {£})) corresponding to the sites 
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m . . . tim, and by u±, . . . vm the rest of the parameters so that {A} U {£} = {[i} U {u}. The 
first matrix element in the formula for Sm is equal to 

({00..0} Jvf {ll..l} M |S(// 1 ) . . . B(ji M )\0) = ({ll..l} M |S , (/i 1 ) . . . BV M )|0), (14) 

where B'(fi) are the same operators difined on the lattice consisting of M sites with the 
inhomogeneity parameters t\ . . . . The second matrix element can be reduced to the same 
expression with the parameters Ai . . . Am (using the symmetry C B, 1). Then using 
the formula for the matrix elements of the operator f^ 1 (see eq.([5p) we finally obtain the 
expression: 

S M ({X},{t})= E Ml a N $ M(t,^M(^)II - f l _ y (15) 
ni,..n M \ j J i,j C v/ i * v 3) 

where we denoted by $m(£, 1) the functions in the right hand side of eq.fllij). The determinant 



representation of this function and its properties are well known (see for example ||,[12 
1,0): 



*) = n (t A n p p\ det (t p m p I \ ^ 

for the rational case. Here {£} are the inhomogeneity parameters and {t} are the arguments 
of B- operators (see eq.([L4])). Note that in the expression fll5|) all the dependence on the 
inhomogeneity parameters £j is contained only in the functions a(i/j). The functions a(^) 
in eq.(|T5|) are exactly the functions defined above a(u) = T\ a c(P, a — v ) while in the rest of 
this formula due to the definition of the matrices T, one should interchange the arguments 
in the functions c _1 (z/j — fj,j) which is taken into account in eq. (|i~5l) (or one could make the 



replacement rj — ► —t]). Using the properties of the functions $m one can represent the general 
formula (|i"5|) in a different way: 

Af 1111 
£ £ (II (a(A„)a(4))) $ m (t fe , A /3 )$ A /_ m (A n , £ a ) - 

m=0 fc,n 



c(A /3 - A n ) c(t Q - 4) c(t Q - A n ) c(X /3 - t k ) 



where the sum is over the two sets hi, ...k m and ...n M _ m . We used the following simpli- 
fied notations in this formula. We devided the set {£} into two subsets {t} = {t k } U {t a } 
where {t k } = (t kl ...t km ) G {i/}, {t Q } G {//} and analogously {A} = {A„} U {A^}, {A n } = 
(A ni ...A„ M m ) G {u}, G {//}. The products in the last formula are over the indices la- 

beling the elements of the corresponing sets. One can further rewrite this formula using the 
expressions for the functions $ m to get the formula which leads to the determinant represen- 
tation for Sm after the special dual fields are introduced [TT|. 



Let us derive the scalar product of the Bethe eigenstate with an arbitrary dual state along 
the lines of ref.0 starting from the formula (|l|). The set of the parameters {t} obey the 
Bethe equations so for each a(U) in the sum (|T5|) one can substitute the function 



ei, \ rr c(t a tj) -j-j- <p(t a tj Tj) 
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Then the idea is that one can calculate the sum in (15) for an arbitrary smooth function a(A), 
which behaves at least as a constant at infinity, used for the terms a(Xj) (not necessarily equal 
to a(A) = Yia c(£a — A)). In that case the sum has the simple poles in the variables {t} and 
{A} at the points tj = Aj and don't have any other poles, for example the poles at Aj = Xj, 
which exist if, according to || one considers a(Xj) as an independent variables. These poles 
was not considered in |§ (actually only the behaviour of the function in the parameters ti 
was considered - obviously there is no poles at tj = tj if the Bethe Ansatz equations are taken 
into account). One can use the symmetry of Sm in {t} and {A} and single out the residue at 
Ai — > ti (which is contained in the term with Ai G {u}, t\ G {/i} and vice versa): 

S M ({A}, {t}, ofa)) | Al _ fl - / (A , l} ~ II * , wx * t\ S "-i « A >'' ' 
where {A}', {£}' do not contain Ai, ti and 

a (Uj) = a{vj)f 



£(ti -Uj)' 

The functions /(ti) entering the sum Sm-i are also modified so that the terms with t\ are 
absent in their definition. Also the behaviour of Sm ~ 1/Ai at Ai going to infinity should be 
taken into account. It can be easily proved that the same recurrence relation at Ai — ► t\ is 
obeyed by the following expression: 

S M ({X},{t}) = T^—r, . * 7T rT det(M 4 ,(t,A)) 

Ui<j{ti - tj) UjKiiXi - Xj) 10 

where at the end of calculations one should take the actual form of the function a(X). So 
we proved that for an arbitrary function a(A) the residues of the only existing poles at 
ti = Xj satisfy the the same recursion relation in both formulas. Since the function (in 
each variable) is completely determined by the positions of (all) poles and the corresponding 



residues (along with the behaviour at the infinity) the two functions ([L7|) and (]T5|) should 
coincide for an arbitrary function a(A). To complete the proof, one should check the formula 
(|TT| ) for M = 1,2. Equivalently, one can check the coefficients corresponding to the terms 
I\iii a (Xi) and ni^io(ti) = 1 in eq. (p7|) which are in agreement with (|15|). The formula (|17|) 
is written for the rational case. This formula can also be represented through the Jacobian 
detij(d / dtiA(Xj] {£<*})• The orthogonality of two different Bethe eigenstates was shown in 
ref.(Q). From the equation (|T7|) taking the limit Aj — > tj one can easily obtain the formula 
for the norm of the Bethe eigenstate M. 



4. Direct proof for the scalar product with Bethe eigenstate. 
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Here we present the direct and the simple proof of the determinant expression ([L7]) for the 
scalar product with one Bethe eigenstate. Although the proof given in the previous section is 
correct the generalization of this method to the expressions for different correlation functions 
is not straightforward. It is still an open question if it is possible to obtain the determinant 
expressions for various correlation functions (the simlest example is the emptiness formation 
probability, which can be represented as a sum very similar to that of eq . (|15|) ) without the 
auxiliary dual fields |f2|,|ll||. Thus the various direct derivations of the formula ([17]) are worth 
of developing. Let us present the direct proof of this formula based on the expression ([15]). 

First, let us rewrite the expression (]T5|) for the case when the parameters {t} satisfy the 
Bethe Ansatz equations. Expressing the functions a(t k ) = f(t k ) and using the explicit form 
of the function $m we get: 

Sm(A, t) = 1 1 E("l) Pfc (-l) P " IIWAn)) 

LLi<j{H - tj) lL,<i(A - A iJ k ,n 

det(A /3 ; t k ) det(t Q ; A n ) \{{Xp - A n + r])(t k -t a + i])(t a - A n + r})(\p - t k + 77), 

where again the products is over all elements of the corresponding sets of parameters. Here 
the notations are exactly the same as in the previous section. The sign factors P k , P n depend 
on the sets of the coordinates k±, . . . k m , n\, . . . nu-m in order to obtain the factors in front of 
the sum from the terms \[{t k — ^a) -1 and II(A/3 — A n ) _1 We also used the simplified notations 
for the determinants: 

det(t, A) = det (- !> ) . (18) 

The order of lines and columns in the determinants det(A,g;tfc), det(t„; A n ) corresponds to 
the prescription k\ < A>2 < • • • < k m and the same for the other sets of variables. Another, 
equivalent way to rewrite the last formula for Sm is: 

= 1 1 n E(-i) ft (-i) f -(-irnw) ( i9 ) 

" An + tj) m -*P-V) det(A,; t k ) det(t a ; A n ) n f **"!"^ ) II i ^'^V ) > 

where again the products is over all elements of the corresponding sets of parameters. Here 
again the sign factors P k , P n depend on the sets of the coordinates k\, . . . k m , n\, . . .nu-m 
and we also used the simplified notations for the determinants (^) and denote for example 
~[\(t — \ n + = Y\j YliL\{U — Kij m this formula the poles at tj = Xj are contained in the 
determinants. If the two last products are dropped out this formula becomes the result of the 
decomposition of the determinant in eq.(|j]). Then one can see explicitly that the residues in 
the poles in the formulas (pl|) and (|17D satisfy the same recurrence relations. 



Our proof of the formula (]T7| ) is based on the on the expression (^). The idea is to 
calculate the sum over the set of the coordinates {k}. For example at the first stage consider 
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the case m = M — 1 which corresponds to the terms containing the single function a(X ni ) (the 
set {X n } contains only one element X ni and the set {t a } - only one element t ai ). Then the 
sum over the set {k} in eq.flTDD has the following form: 

.(A., ) g dett V U) K _ _ _ n (j^±g) n . (20) 

This expression can be considered as a first column developement of the determinant of some 
new M x M matrix, the parameter A m being assigned to the first column of this matrix. 
To find this matrix let us first consider the determinant 

/ i \ 

det (Mij) = det 



(U - Xj)(U - Xj + T])J ' 

and modify the first column of the matrix Ma adding a linear combination of the other 
columns so that the determinant det M is unchanged: 

Let us choose the coefficients as: 

fei,x ~~ V) a (Ai -t a -rj) 

Notice that this expression does not depend on the index i. It is necessary to calculate the 
following sum: 

where 

f M \ _ TT (*a - X x + 77) 

AArJ - 11 77 : — ; — r- 

To calculate the sum in eq.(^) consider the integral in the complex plane over the circle of 
large radius which is equal to zero due to the behaviour of the integrand at infinity: 

dz M = f(z) = TT ^ a ~ z + r]) 

The sum of the residues at z = Xp is equal to the sum in (pl|) , so it is easily calculated to be 

/(Ai) + f(ti) 



(Ai - fc) LVi(Ai - A^) (t, - AO n^i ft - A/?) ' 

Thus combining all terms we obtain the following expression for the new elements of the first 
column: 

M f = 1 TT ( A l - V) T-r (ta -tj+rj) 
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Note that there are several ways to choose the coefficients C x in order to obtain the closed 
expressions of this type for the first line or column for the matrix M. 

Now we are ready to come back to the sum (PU[). Let us define \'p = \p + rj and take into 
account the equation det(A / 3; tk) = detf^; A^). Then we get the following sum: 

a( A„,) B-ir dette 4)7— ^——n ( 7^4^) n ( x 4~ M 



ai (tai ~ X ni )(t ai — A ni + rj) \tk — A ni + r] J y \'p — t ai J 

Comparing this sum with the last formula we find that this sum it equal to the determinant 

a(A ni )det({t};{A„ 1 A^}), 

where the parameter A ni corresponds to the first column of the new matrix. This determinant 
in turn can be represented as its first column developement: 

«(A ni ) E(-i) ai 7T 3— — 1 — r~^ det (A/3; tk). 

ai [t ai — Arnica! — \ n -T V) 

Comparing this expression with (^) we see that the two last products in (^) can be removed 
- the sum over a>\ does not change. Now let us come back to the formula (|T§D. Instead of the 
the summation over the set {t^} one can take the sum over the set {t a }. One can label this 
set by the coordinates aii, . . . au- m on the lattice consisting of M sites with the corresponding 
parameters t\,t% % . . Am- We can consider a« as independent coordinates on the lattice, the 
only restriction being oti 7^ ctj. The determinant det(t a ; A n ) can be represented as the sum 
over the permutations of the products of matrix elements corresponding to the pairs of the 
variables t a . , A Wj . Each term in the sum contains all the variables t ai and all the variables X nj . 
For each term one can take the sum consequently over each of the coordinates ai, . . . otM-m 
with all the other indices fixed, so that only the set {k} changes during the summation at 
each step of this procedure during the summation over each ai. For each ai one can single out 
the terms with ai and the corresponding rij from the last two products in the formula ([19|). 
Then we perform the procedure described above. One can see that for each term in the sum 
over the permutations (for the determinant det(t Q ; A n )) this procedure removes consequently 
all terms in the last two products in the formula (0) and we obtain: 

lli<j\ l i l j) llj<i\ A i A j) k,n 

H(t - A n + ri) Y[(t -Xp-ij) det(A /3 ; t k ) det(t a , X n ). (22) 

This formula is the result of the decomposition of the determinant of the sum of the two 
matrices in (|T7|). In fact, due to the symmetry of Sm{X, t) in the variables Aj, it is sufficient 
to prove this statement for the set {no} = (1, . . .p),p = M — m, and {/?o} — (p + 1> • • ■ M). 
such that (— l) Pn o = 1. Clearly, the factor (— l) Pn makes the expression for Sm symmetric 
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with respect to {A}. For the term with {n} = {no} the determinant in eq.fllTD is represented 
as a sum over the permutations P G Sm with the sign (— l) p : 

S M {\t)\ {no} = - — — — — T-rl[a{\ no )f + {\ no )l[f-{\p ) 

(_l)m J2(-l) p J] a(t Pno , A„ ) n a(\ Po ,t PPo ), (23) 
p 

where / ± (A J ) = Ila(^a ~~ A j ^v) an d a (^; A j) are the matrix elements (|I~8|). Each permutation 
P can be represented as P = P(a)P(k)P({n }, {«}) where the permutation P({n }, {a}) is 

(1,2, . . .p,p + 1, . . . M) — * (ai,a 2 , . . . a p , . . . k m ) 

(here it is implied a, < a i+ i, fcj < fcj+i) and the permutations -P(a), -P(fc) corresponding to the 
permutations of indices a, and fc, among themselves produce the determinants det(t a ; A no ), 
and det(A ( g ; tk). The sign of the permutation P({n },{a}) differs from the sign (— l) Pfe by 
(— l) pm since (— l) Pfc corresponds to the permutation (fci, . . . k m , a±, . . . a p ). Thus comparing 
the signs in the equations (|23|) and (p2[) we see their equivalence. Thus we reproduce exactly 
the formula (|17|). 
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